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ABSTRACT 

We describe a microfluidic device for studying the orientational dynamics of microrods. The device enables us to 
experimentally investigate the tumbling of microrods immersed in the shear flow in a microfluidic channel with 
a depth of 400 (i.m and a width of 2.5 mm. The orientational dynamics was recorded using a 20X microscopic 
objective and a CCD camera. The microrods were produced by shearing microdroplets of photocurable epoxy 
resin. We show different examples of empirically observed tumbling. On the one hand we find that short stretches 
of the experimentally determined time series are well described by fits to solutions of Jeffery's approximate 
equation of motion [Jeffery, Proc. R. Soc. London. 102 (1922), 161-179]. On the other hand we find that 
the empirically observed trajectories drift between different solutions of Jeffery's equation. We discuss possible 
causes of this orbit drift. 
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1. INTRODUCTION 

In recent years there has been a significant development of lab-on-chip devices designed to be used as biological 
assays^ in environmental monitoring for the detection of toxic materials^ as well as for the assembly of nano- 
and microscale objects into more complex systemsEE] The main reason for miniaturisation is to reduce the cost 
of production for the devices, as well as to reduce the volume needed to analyse samples. Last but not least, 
the lithographic method used to make the devices allows the production of arbitrary channel structures at a low 
cost. 

There have been numerous investigations where microfludic devices have been used to study the properties 
of microparticles. For example, Eriksson et al. used a microfluidic system in combination with optical tweezers 
to investigate salt stress in yeast cells!^ Alrifaiy et al. conducted electrophysiological investigations of individual 
cellsP Tegenfeldt and coworkers have developed a method to sort cells that differ in both size^and morphology 8 
by means of their motion in an array of obstacles. These studies rely on a precise understanding of the dynamics 
of particles in microfluidic devices. The translational dynamics of microparticles is often well described by 
center-of-mass advection. 

For non-spherical particles the orientational dynamics must also be considered. Already in 1922 Jeffery de- 
rived an approximate equation of motion for the orientational dynamics of non-spherical particlesP Jeffery's 
solution predicts that elongated axisymmetric particles mostly align with the flow direction, but periodically 
rapidly turn around 180 degrees. This motion is referred to as 'tumbling' in the literature. Jeffery's equations 
form the foundation for studying the rheology of suspensions of non-spherical particles, together with a diffusion 
approximationEnHH] We refer to Petrie^l for a contemporary review of fibre-suspension rheology and the appli- 
cations of Jeffery's equations in this field. Recently, pattern formation by elongated particles in random flows 
was investigated using Jeffery's theorypKIl with applications to flow visualisation (rheoscopy). With much past, 
present and future research based on Jeffery's equations, it is necessary to study how accurate Jeffery's theory 
is in predicting the orientational dynamics of a single microparticle. 

The aim of this investigation is to design a device that makes it possible to study the orientational dynamics 
of a single non-spherical particle, such that the results can be compared to theoretical predictions. In this article 
we describe a microfluidic device allowing us to study single-particle orientational dynamics advected in a channel 



flow. The channel is produced in PDMS, and the particles are produced using a liquid-liquid dispersion technique. 
By video microscopy with a motorised stage we observe the orientational motion of particles advected along the 
channel. A recent experiment)^ displays short sequences of tumbling motion of Escherichia coli bacteria. It 
is demonstrated that the bacteria follow Jeffery's theory as they tumble a few times. Our experimental setup 
enables us to follow a particle over longer distances. We show examples where stretches of the experimentally 
determined time series are well described by fits to solutions of Jeffery's equation of motion. But we also 
demonstrate that the empirically observed trajectories drift between different solutions of Jeffery's equation. We 
discuss possible causes of this effect: thermal noise, random fluctuations of the flow due to imperfections in the 
channel, and the possibility that the particles are not perfectly axisymmetric. We argue that inertial effects, due 
to particle or fluid inertia, are less important. 

We see the work described in the following as a step towards a device that makes it possible to investigate the 
dynamics of a single particle with full parameter control. Our aim in future studies is to control the particle shape, 
its initial conditions, and the flow conditions. Particles produced with lithographic 17 and photo-polymerisatiorP^ 
methods are ideally suited to achieve precisely tailored particles. The initial conditions of a particle, both with 
respect to position and orientation, can be controlled by means of optical tweezers. The viscosity of the medium 
is conveniently varied by using different concentrations of a glycerol/ water mixture. 

This article is organised as follows. In Sec. [2] we describe the experimental setup, the fabrication method for 
rod-like microparticles and the data analysis procedure. Our measured experimental data is presented in Sec. [3] 
The results are compared to theory and discussed in Sec. |4j We end the article with concluding remarks and an 
outlook in Sec. [5] 

2. MATERIALS AND METHODS 

2.1 Microrod fabrication 

The polymer microrods used in this experiment are produced by a liquid-liquid dispersion technique.^ A 
photocurable epoxy resin, SU-8, diluted in solution of 7-butyrolactone (GBL), is introduced into a 1:1 mixture 
of ethylene glycol and glycerol while stirring at 600 rpm. Using a syringe, the SU-8 solution is injected into 
the spinning liquid where it forms microdroplets. These droplets are quickly drawn out to urn-sized rods by 
the shear forces in the spinning glycol/glycerole mixture. As the GBL transfers from the SU-8 solution into 
the surrounding mixture, the resin solidifies in the shape of microrods. After 10 minutes of stirring at constant 
rotation the microrod solution is introduced to irradiation of ultraviolet radiation with a wavelength of 365 nm 
for 15 minutes to experience an irreversible cross linking of the photocurable resin. Thereafter, the polymer 
microrods are separated from the non-solidified polymer solution using repeated dilution with deionised water 
and centrifugation. Finally, the microrods are suspended in a solution consisting of a mixture of water and 
glycerol. The polymer microrods have a tendency to entangle or coalesce in the solution. This problem is solved 
by adding a few percent of surfactant to the solutions, which reduces the effect of entanglement and coalescence 
Fig. [IJl shows a transmission image of the produced microrods using a 20X microscope. The rods produced using 
the method described above have a high aspect ratio with a rather wide size distribution. Fig. [TJ) shows the 
measured size distribution in a sample of microrods. 

2.2 Design and fabrication of microfiuidic channels 

A single-inlet /single-outlet microfiuidic device (2.5 mm wide and 400 |i.m deep) is produced in PDMS. A molding 
form is first created in aluminum using a high precision milling machine. PDMS is molded in this form, cured, 
removed and finally hole-punched. The PDMS is irreversibly sealed to a plasma treated glass microscope slide 
with a thickness of 1 mm. This relatively thick glass is selected over the standard cover glass in order to avoid 
mechanical deformation of the cover glass caused by the PDMS. Initial work using a thin cover glass showed that 
the cover glass deformed. This caused the images of the moving particle to get out of focus as it was tracked 
along the microchannel. 
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(a) A transmission image of polymer mi- (b) Measured length distribution of parti- 

crorods imaged through a 20X micro- cles. The particles used in this work are ap- 

scope objective. The distance between proximately 30 [im long, 
minor ticks is 10 [im. 



Figure 1: Polymer microrods. 
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Figure 2: The experimental setup with coordinate system. 
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Figure 3: Three frames from video microscopy. The images are in sequence and the time between subsequent 
frames are approximately one second, corresponding to 100 frames. The image is 348 x 260 pixels wide and high, 
respectively. Microrods move from left to right, and the camera is following. The particular microrod in focus is 
indicated by an arrow in each frame. 

2.3 Experimental setup and procedure 

Our experimental setup is schematically shown in Fig. [2] Microrods moving in the microfluidic channel are 
manually tracked using a motorised stage. The viscosity of the fluid is selected by using a mixture of water and 
glycerol. A 1:1 ratio is used in all the experiments described below. The water /glycerol mixture is pumped into 
the microchannel with a calibrated syringe pump (CMA 400/ Microdialysis) at flows ranging from 5 ul/min up 
to 11 |j.l/min. A 20X microscope objective (NA = 0.28) is used to image and record the microrods' motion onto 
a CCD camera (Leica DFC 350 FX). The acquired videos are analysed using MATLAB (Mathworks, Natick, 
MA, USA). 

2.4 Data analysis 

The camera used in the experiment produces a gray-scale video captured at 100 Hz at a resolution of 348 x 
260 pixels. A typical sequence of one particle during its movement in the channel is several minutes long. Fig. [3] 
shows three frames from one experiment presented in this article. The main observable is the orientation of the 
particle as a function of time. The orientation of an axisymmetric particle can be described by a single vector 
along the particle's symmetry axis. We define the unit vector n(i) as the direction of the particle's symmetry 
axis at time t. 

The particle orientation is described in a Cartesian coordinate system where the x-axis is along the channel, 
the y-axis pointing down into the channel, into the camera, and the z-axis perpendicular to the channel direction 
in the camera plane. The coordinate system is shown in Fig. [2] It is important to keep in mind that the video 
captures a two-dimensional projection of the three-dimensional system. The component of n(t) in the y-direction 
(see Fig. [2]), can therefore only be determined in magnitude, but not in sign. 

We employ the following procedure to analyse the video data in two steps. 

1. The endpoints of the particle are identified in each frame. 

2. The endpoints in the camera plane are used to compute the orientation vector n(t). 

The first step, identification, is performed by assuming that we know the endpoints in the previous video frame. 
A neighborhood of the particle is cropped from the video and converted to a binary, black-and-white image 
by Otsu's method PS All remaining connected domains are fitted to ellipses by means of principal-component 
analysis, and the ellipse closest to the previous frame is chosen. Then the video is advanced one step, whereafter 
the procedure is repeated. This procedure may fail for several reasons. For example there may be no connected 
domains in the binary image, or the ellipse closest to the previously found may in fact be the wrong one. On 
failure the program stops the automatic identification. In this case the endpoints of the rod must be determined 
manually and entered into the program. 
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Figure 4: Components of the unit vector n(t) as a function of time. The top panel shows the x-component, which 
is along the flow in the channel. The middle panel shows the magnitude of y-component, which is perpendicular 
to the camera plane. The bottom panel, finally, shows the z-component, or spanwise in the flow. The value of 
n z as the microrod passes through zero yields a measurement of which particular orbit the microrod is following. 

The principal-component analysis cannot distinguish the difference between the two ends of the particle. It 
only gives the angle of inclination in the interval [—90, 90] degrees. Therefore we manually mark where the 
direction defined by the endpoints actually pass the 90 degree inclination. After this step we have acquired a 
time series of the distinguishable endpoints of the particle. 

The particle projection of n(t) onto the camera plane and the length of the microrod can then be used to 
determine the full orientation vector n(t), up to the sign in the y-component which remains undetermined. 



3. RESULTS 

The experimental data obtained in this work are time series of the unit vector n(i) describing the microrod 
orientation. An example is shown in Fig. [4] The x-, y- and z-components of n(i), denoted n Xl n y and n z , are 
shown separately. The x-direction is along the flow and typically stays close to n x = ±1 with periodic switches in 
sign. This type of motion is referred to as tumbling motion, as explained in the introduction. The y-component 
is perpendicular to the camera (x — z) plane. It has to be measured by variations in the length of the microrod 
in the camera plane which presents two problems. First, it only gives information on the magnitude of the 
y-component. Second, the signal is very noisy when the microrod is aligned in the x — z-plane, which is the case 
in a large part of the time series observed. 



4. DISCUSSION 

In Fig.[| we observe long periods of alignment of the orientation vector n(i) along the flow direction, interspersed 
regularly by tumbling events where the particle quickly turns 180 degrees and aligns in the opposite direction. 
Between such events we have \n x \ « 1, and consequently both n y and n z are close to zero during these periods. 
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Figure 5: Schematic drawing of the flow profile u(y). Note that the view shown in this figure is perpendicular to 
the view seen through the microscope objective. The left part of the figure shows the flow profile over the depth 
of the channel, with the tangent if the profile indicated at the position of a particle. The right part of the figure 
shows a linear approximation of the flow in a coordinate frame translating with the particle along the channel. 



During the rapid tumbling phases, typically both n y and n z vary. But we have also observed motion purely in 
n y . 

The following paragraphs describe the quantitative comparison of our results to Jeffery's theory. Jeffery's 
results rest on the calculation of the hydrodynamic forces and torques on an ellipsoid, in the limit of a viscous 
linear flow. We use his results for a special case with three assumptions. First, the ellipsoid is axisymmetric. 
Second, the particle's inertia is negligible. Third, the flow is a shear flow, as depicted in Fig. [5] For this case. 
Jcffery showed that n(t) must describe one of infinitely many closed orbits on the unit sphere. Which particular 
orbit a particle follows is determined by the initial condition. The complete dynamics is determined by the shear 
direction, the shear strength s, and by the aspect ratio A of the microrod. Examples of orientational dynamics 
are shown in Fig. [6] When the particles are not perfect ellipsoids, A must be interpreted as an effective aspect 



ratio (see the discussion in Ref. 22 1. 



The channel used in this experiment has an aspect ratio of 1:6.25. We can therefore assume that far away 
from the side walls, there is no component of the flow gradient in the z-direction. Further, the particles are short 
(order of 10 (a.m) in comparison to the channel depth (several 100 urn). Thus, locally the particles experience a 
shear flow in the x-direction, with shear in the y-direction and vorticity in the z-direction. Fig. [5] depicts the 
linearisation around the particle. In this approximation the flow field u is locally described by its gradient 




A = Vu= , (1) 



where s is the shear strength. Jeffery's equation for an axisymmetric particle is 

^=Bn-n(n T Bn), with B = — ^—-r(X 2 A — A T ). (2) 
at 1 + \ l 

It can be showrEHEIl that this non-linear equation is solved for any time-independent matrix A by 

q(i) = exp(flt)n(0) and n(t) = , (3) 

where n(0) is the initial orientation at t = 0. If we assume that n x — at t = and insert the flow gradient (JTJ 
into the solution ([3| we arrive at 



. I A sin uit\ . 

n(i) = — cos ujt , uj = -. (4) 

Vtan 2 6» + cos 2 uit + A 2 sin 2 cot \ tan g Q ' 1 + - v 
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(a) Jeffery orbits 



(b) As seen from camera 



Figure 6: Visualisation of Jeffery orbits from Eq. Q with A = 10 and 6 = ^Xi, Xi = 0, 0.3, 0.6 from bottom to 
top. (a) illustrates the closed paths n(t) takes on the sphere where the dots are plotted for equidistant times, (b) 
shows the same but from the viewpoint of the camera in our experiment. A measurement of n z at the midpoint, 
where n x — 0, uniquely determines the orbit of the particle. 

Here 9o is the angle between n(0) and the z-axis when n x — 0, effectively implementing an initial condition and 
choosing an orbit. The angle 9q is therefore referred to as the orbit constant^ Examples of the solutions |4]) are 
shown in Fig. [6] 

In general any point n on the unit sphere is associated with a unique orbit. In our particular case the most 
convenient point to measure which orbit a particle is following is at n x = 0, in other words in the middle of a 
tumbling event. Reading off the peak height in n z where n x passes through zero, as indicated in Fig. |6b| uniquely 
determines the orbit the microrod is following. 

Using the z-component of Eq. Q we find segments in our data closely following a Jeffery orbit. We show 
time series of n z with such excerpts in Figs. [7j and [5J Fig. [7] depicts a sample where the particle stays on a similar 
trajectory for several consecutive tumbling events. The sample in Fig.[HJ by contrast, shows deviations which we 
interpret as orbit drift between subsequent tumbling events, as further discussed below. 

There are three free parameters to consider when Eq. Q is used to describe the motion of the particle. They 
are the shear strength s, the effective aspect ratio A, and the orbit constant 8q. The following paragraphs discuss 
these parameters in order. 

The shear strength s affects the period time of tumbling linearly. The shear strength s is determined by the 
flow profile in the channel, and by the position of the microrod in the y-direction. Moreover, the shear rate can 
fluctuate if the syringe pump does not provide perfectly steady flow. Currently temporal variations of the shear 
rate are not available with good enough precision, and s must remain a free parameter. As shown in Figs. [7j 
and [5J the estimated s is of order unity. This is consistent with the setup of our experiment where the depth 
of the channel is L s» 400 u\m, making the typical flow velocity v ps sL/2 = 200 [iiri/s. The channel is around 
5 cm long, and it takes the particles around 5 minutes to travel its entire length. A speed of 1 cm/minute gives 
approximately 150 (im/s. The observed speed is hence consistent with the estimated shear strength. 

The effective aspect ratio A is a dimensionless number describing the shape of an axisymmetric particle. An 
aspect ratio of A = 1 corresponds to a sphere, and A = 10 corresponds to a cigar-shaped particle. For an ellipsoid, 
A is the actual aspect ratio. For a microrod we expect A to be close to its length divided by its width. The 
estimates in Figs, [jj and [8] are A = 6 and 7, respectively. For both examples presented we measure A ~ 10 (see 
Fig. [3]), in reasonable agreement. The aspect ratio is a particle property and as such it cannot be time dependent. 

*The common notation for the orbit constant in the literature is C (see for example Eq. (3b) in Ref. |24[) which relates 
to our notation by C = cot So- However, we prefer #0 as it is a bounded variable on [— f , f] with a clear physical 
interpretation as an angle. 
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Figure 7: Sample trajectory showing a relatively stable orbit over time, (top) Component n z (t) extracted from 
a video sequence displaying numerous flip events on very similar paths, (bottom) Excerpts from top image, 
showing in detail the comparison to a suitable Jeffery orbit (dashed) from Eq. Q. The parameters for the fit 
are A = 6 and for the three excerpts s — 1.80, 2.00, 1.95, respectively, and orbit constant 6q — 0.49, 0.46, 0.51. 



The orbit constant do is a free parameter, a constant of integration, in the derivation of Eq. Q. The theory 
gives no preference to any particular orbit, but asserts that O must be constant in time. In the figures, and 
as can be seen in Eq. Q, 9$ determines the peak height of n z as discussed above. Clearly the observed values 
of 6q are not perfectly independent of time. In the context of Jeffery's theory this may be described as drift 
between different solutions of Eq. Q , or orbit drift. We infer that there must either be effects beyond the purely 
hydrodynamic forces considered by Jeffery (possible effects may be inertial, due to either particle inertia - finite 
Stokes number - or fluid inertia - finite Reynolds number). Another reason may be that the particles are not 
perfectly axisymmetric. Last but not least the particle dynamics is certainly affected by noise, either thermal 
noise or due to flow irregularities. These possibilities are discussed in the following paragraphs. 

A fluid Reynolds number based on the channel dimension (order of 10 2 (J.m) and the flow speed (order of 
10 2 \im/s) is on the order of Re = 10 2 , for water at room temperature with kinematic viscosity of 10 6 m 2 /s. 
In our experiments the viscosity is even higher since we mix glycerol into the liquid. The particles used in this 
work are neutrally buoyant. The Re-dependence has been studied numerically by several authors JHEH Their 
results indicate that the Jeffery solution based on creeping flow (Re= 0) is valid well up to Re = 10 _1 and 
possibly higher. Lundell and Carlsson^l studied numerically the effect of increasing St at low Re. They found 
that the particles follow the Jeffery dynamics for moderately small St with only a very slow orbit drift towards a 
neutral center orbit (corresponding to 9 = in our notation). On the contrary, we observe that orbit drift can 
occur quite rapidly, and in any direction, not towards a particular equilibrium. Hence, previous investigations 
of inertial effects do not provide an explanation for our empirical results. 

The effects of noise on the single-particle orientational dynamics have not received much attention in the past. 
However, recent years have seen an increased interest in single-particle dynamics, sparked by single-polymer 
experiments such as those by Smith et aiP^ and Schroeder et al., 29 followed by theoretical investigations by 
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Figure 8: Sample trajectory showing clear orbit drift, (top) Component n z (t) extracted from a video sequence 
displaying numerous flip events with varying paths, (bottom) Excerpts from top image, showing in detail the 
comparison to a suitable Jeffery orbit (dashed) from Eq. Q. The parameters for the fit are A = 7 and for the 
three excerpts s = 3.40,4.60,3.80, respectively, and orbit constant 9q — 0.57,-0.85,0.57. 

Chertkov et alW^ and TuritsynP^ These studies focus on the limit of extremely slender particle^ A — > oo, 
and assume that there is a random component of the flow gradient, A, modeled by a diffusion coefficient. As it 
turns out the finite aspect ratio A is crucial to the Jeffery dynamics of periodic tumbling. Since our observations 
suggest that the tumbling motion is indeed periodic, the limit A — > oo is irrelevant to our experiment. 

There have been several studies of angular statistics of large ensembles of particlesJ-SED However, these 
studies focus mainly on the rheology of macroscopic suspensions. We can use the expression in Ref. [TT] to 
approximate the rotational diffusion constant of a spheroid, assuming room temperature and a 10 x 1 x 1 urn 
particle in water. We find that D r s» 10~ 2 radians/s. The rotational diffusion rate should be compared to 
the typical rate of rotation due to the shear. While the particle is nearly aligned with the flow direction the 
deterministic rate of rotation is of order s/(l + A 2 ). With s»l and A « 10 we find that the deterministic and 
noise-induced rates of rotation are comparable in magnitude. This discussion leads to the conclusion that noise 
is an important mechanism to describe our observations. 

Finally, we consider the possibility that the particles are not perfectly axisymmetric. Even small deviations 
from axisymmetry can result in relatively large changes in a particle's rotationEHEII In particular the results 
of Refs. [32] and [33] indicate that weak asymmetry can result in drift between different observed values of the 
orbit constant. In order to experimentally distinguish between the effects of noise and particle asymmetry, it 
is necessary to precisely control the particle shape. Furthermore, observing more consecutive tumblings will 
facilitate to distinguish periodic modulations due to asymmetry from random drift. 

5. CONCLUSION AND OUTLOOK 

To conclude, our results suggest that a plausible way to understand the empirically observed orientational 
dynamics of single particles may begin in the hydrodynamic theory by Jeffery. We have argued that two additional 

^ Actually, the authors investigate flexible polymers, which do not have the same dependence on the aspect ratio. The 
equations, however, are equivalent to the limit of rigid rods of infinite aspect ratio. 
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mechanisms may affect the dynamics significantly. First, noise causes random drift between different Jeffery 
orbits. This noise could be modeled as a random contribution to the flow-gradient matrix A. The fundamental 
parameters to explore in this case are the aspect ratio A, and the Peclet number. The Peclet number, as is well 
known, is a dimensionless ratio between the magnitudes of the deterministic and random components of the flow. 
The second mechanism that deserves further investigation is that weak asymmetry of particles can induce drift 
in the orbit constant. The case of non-axisymmetric particles is part of the Jeffery theory, of which Eq. ^ is a 
special case. 

Future experiments will focus on an increased control of the relevant parameters. Work in progress includes 
coupling a piezoelectric clement to the glass cover, which can be used to induce ultrasonic noise in the liquid. 
With particle velocimetry and volumetric flow measurements we can determine the shear strength. The initial 
position and orientation of the microrod can be controlled by installing optical tweezers in the setup. In parallel 
we are developing a photopolymerisation technique to consistently produce high precision particles of arbitrary 
shape. As we have argued, precise control over the particle shape is crucial to distinguish the effects of particle 
shape from random orbit drift. 
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